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Abstract. We examine the impact of discrete numbers of stars in stellar populations on the results of Chemical 
Evolution Models. We explore the resulting dispersion in the true yields and their possible relation with the 
dispersion in observational data based on a Simple Closed-Box model. 

In this framework we find that the dispersion is larger for the less evolved or low abundance regions. Thus, the 
age-metallicity relation may be a tracer of the Star Formation History of our Galaxy. This theoretical dispersion 
is especially high for the relative abundance log(N/0) in regions where the total number of stars created is still 
low. This may explain part of the scatter in the N/O ratio observed in star forming galaxies. 
We have also found a first order theoretical estimation for the goodness of a linear fit of the Helium abundance 
vs. 12 + log (0/H) with values of the regression coefficient between 0.9 and 0.7 (independent of sampling efi^ects). 
We conclude that it is necessary to include these sampling effects in a more realistic Chemical Evolution Model in 
order that such a model reproduces, at the same time, the mean value and the dispersion of observed abundances. 
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1. Introduction and motivation 

The modelling of any observable related with the inte- 
grated contribution of stellar populations needs to assume 
an Initial Mass Function (IMF) and a Star Formation 
History (SFH). The representation of such distributions 
in the form of analytical expresions allows an easy im- 
plementation in any evolutionary model. However, such 
analytical laws are valid only under the assumption of an 
infinite number of stars, when sampling effects (i.e. the 
integer nature of the number of stars) are not relevant. 
Even in the case of the modelling of a system where the 
real physical properties are known, sampling effects must 
be present showing a dispersion of the data around the 



mean value obtained by the model (see, e.g. Cervino et al 
200C| ). 

Cervifio et al. ( ^002 ), based on the work by Buzzoni 
( |1989| ), have presented a method to evaluate quantita- 
tively the intrinsic dispersion due to sampling effects in 
spectrophotometric synthesis models. It is important to 
remark that the dispersion associated to such sampling ef- 
fects is indeed observed, and used as a distance indicator, 
in the form of surface brightness fluctuations in elliptical 



galaxies (see Buzzoni 1993 , for a theoretical study of the 
subject). 

In a next step, [Cervino et al. ( 2001 ) applied this for- 
malism to some quantities which are time-integrated to be 
obtained, such as the ejected masses of elements by the 
stars when they evolve. With that work, it was demon- 
strated that the estimation of these ejection rates may 
have a large dispersion, mostly in the first 5 Myr of the 
evolution of a stellar cluster. 

On the other hand, the observational data used in 
the comparison of Galactic Chemical Evolution Models 
(GCE) have a very large dispersion: in the Milky 
Way Galaxy (MWG) there exist regions at a same 
Galactocentric distance with oxygen abundances differing 
by 0.2 dex, means dispersions as large as 50%. Such dis- 
persion can be due to observational errors, but also to 
sampling effects. As an example, a stellar cluster that has 
transformed gas into stars following a Salpeter IMF in the 
mass range 120 - 0.08 Mq in a single burst, will produce 
6.8 X 10"'^ SN Mg^. In the galactic context, the most mas- 
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sive O B association known (Cygnus 0B2, c.f. Knodlsedei 
has transformed into stars 6.5 x 10^ M© of g as, as- 
suming the presence of 120 O stars in the association (that 
is, in fact, an upper limit). A region of such characteristics 
will produce 440±20 SN in a 63.8% (1 a) confidence inter- 
val due to the limited number of stars in the region. The 
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formal relative dispersion is around 5%, which will result 
in a similar relative dispersion in the observed mctallicity 
of that type of regions. So, if the chemical evolution of our 
Galaxy is driven by such a type of massive stellar clusters, 
at least, a theoretical relative dispersion around 5% must 
be expected in the comparison of GCE models with ob- 
served data. But in fact, the dispersion must be larger as 
long as most of the observed OB associations have a lower 
amount of gas transformed into stars and the enrichment 
depends on the initial mass of the exploding star. 

Taking into account the dispersion found for the stel- 
lar yields, we think that the same sampling effect may 
also produce a dispersion in the interstellar medium abun- 
dances. However, the stellar yields do not give directly the 
abundances in a region. It is necessary to use them as the 
input in a chemical evolution model. Chemical evolution 
models are usually computed through numerical methods 
which solve a system of equations where the IMF such as 
the SFH are taken into account. However, before comput- 
ing the dispersion of the abundances predicted in these 
numerical chemical evolution models, we would like to es- 
timate the importance of these sampling effects with sim- 
pler chemical evolution models, and see if the resulting 
dispersion is at least of a similar order of magnitude to 
that observed. Therefore, the objective of this paper is 
to obtain a first order estimation of how relevant these 
sampling effects can be in simple GCE models in order to 
consider whether it is necessary to include them in more 
complicated numerical chemical evolution models. 

The structure of the paper is the following: In Sect. 
2 we show the application of the formalism presented in 



Cerviiio et al. (2002) to the computation of yields in GCE 



models. In Sect. 3 we obtain the dispersion on the true 
yields for several metallicities and turn-off masses. In Sect. 
4 we apply the formalism to a GCE closed model and 
we obtain a first order approximation of the relevance of 
sampling effects in the results. We show the conclusions 
and implications of this work in 5. 

2. The computation of stellar ejection rates 
including sampling effects 

Let us assume a system where, at a given time t, there are 
n^{mi, t) — n^^i stars of a given initial mass rrii which have 



completely evolved and nii{mj,t) 



stars of a given 



initial mass rrij which are still present in the system. The 
total number of completely evolved stars, Nj°*-{t), and of 
stars still present in the cluster, Nj°^{t), are: 



Nr(t)^yn^4t) 



Nl°\t) 



i=l 



E 

i=l 



n.At), (1) 



where Id{t) and /a(i) are the number of initial masses 
considered and nd,i(i) and na,i must be integer numbers. 

Let us assume that each star of mass rrii is able to 
eject an amount of freshly synthesized element Z given by 
myz^mi) — myi,i and defined as: 



"T-2/Z,i = "l(Z)ojoctcd,i - z{mi - mi-cm,i) 



(2) 



where m(Z)ojoctcd,i is the amount of the element Z ejected, 
z is the mass fraction of element Z in the star at the begin- 
ning of its evolution and mrom.i is the mass of the stellar 
remnant of a star with initial mass rrii at the end of its 
evolution. 

The amount of freshly mass of the element Z ejected 
into the interstellar medium at the time t, Yz{t) is the sum 
over all the evolved stars of the individual contribution of 
each star: 



id{t) 

Yz{t) = ^ nd^i{t)myz,i 
1=1 



(3) 



The corresponding variance of Yz{t), cr^(Yz(t)), is the 
sum of the individual variances. Applying Poisson statis- 
tics (cr^(nd,i) = rid,i) we obtain: 



Id Id 
o-^(^z) = ^ o-^(»T-d,i) {rnyz^if = ^ nd,i {niyz,i)'^. (4) 



i=l 



i=l 



The amount of material locked into stars and stellar 
remnants is: 



id{t) 



Mt) 



Mock(i) = nd,i{t) TOrom.i + "a,i(i) nij , (5) 



1=1 



The corresponding variance cr^(Miock(i)) is: 
id{t) Mt) 

CT^(Mock(t)) = Yl "-d,z(0 TOrem,* + E "aj (0 (6) 



The yield of the element Z, Pz{t), is defined as the 
fraction of the amount of ejected material over the amount 
of material locked into stars and stellar remnants: 



Pz{t) 



Yz{t) 



Mlock(t) 

It has an associated variance: 



2 _ a^iVz) , r|a2(Afiock) ^ Yz 



<ck 



<ck 



(7) 



cov(Afiock,i"z),(8) 



where cov(Afiock, ^z) is the covariance: 



cov(Miock,>z) 



Id 

E 

1=1 



(9) 



So, at this moment we have the basic ingredients for a 
GCE model and their corresponding variances. The IMF, 
$(m) and the SFH, ^'(i), are implicitly assumed in the 
computation of the ni{t) values and so, the yield. Since our 
objective is obtain a first order estimation of the relevance 
of the discrete number of stars (i.e. sampling effects) , we 
will only examine the case of the true yields. 
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3. True yields 

The true yield is defined as the resuhing yield of a sin- 
gle stellar population. It is in fact, the yield obtained 
from a generation of stars that have suffered a single 
Instantaneous Burst of star formation: '^{t) — Ms x S{t) 
where Alg is the amount of gas that has been transformed 
into stars in the IMF mass range m^p and miow In this 
situation the values of ni{t) can be obtained by a direct 
integration of the IMF: 

/■"Kip /-"lup 

Nj°\t) ^ Ms / 4>(m) dm = Ms A m" dm (10) 

{t) <y (t) 

where m^p is the upper limit of the IMF, mr{t) is the 
initial mass of the star that ends its evolution at time 
t, that can be approximated by a power lawQ: mr{t) = 
Bt~^. The IMF has been normalized to 1 M© with limits 
TTT-up and miow with a normalization constant A. It has 
been assumed a power-law with slope a over all the mass 
range. Using this approximation, Nj°^ has the functional 
form: 



N, 



Tot 



(11) 



Let us assume that myz.z has a functional form 
myz{m) = am + b. In this case, Yz{t) and a'^{Yz{t)) will 
have the functional form: 



Yztt) - Ms 



cp(m)(aTO H- 6) dm 



M,Ai 



-,a+2 



a + Z 
Ms A b 
a -I- 1 



(12) 



Ms / ^{m){am + b f dTTi 
M,Aa^ 



a + 3 

MsA2ab . .R,-7^a+2^ 
M Ab"^ 



(13) 



In the same way, if we assume that m-^c-a^[m) has also 
a functional form: mrcm(TO) = cm + d, we may compute 
easily Miock and (T^(Miock)- Both turn out to be functions 
depending on m^p, miow and t. 

Let us now show how the relative dispersion in the dif- 
ferent quantities scales. At this point it is useful to define 
it in terms of J\f = A^off- This quantity, already used in 
previous works of this series (see also 
defined, for any quantity A as: 



Buzzoni 



1989), is 



A 



1 



(14) 



^ However, see [Cervino et al. ( 2001 ) for remarks in this ap- 
proximation. 



It means that the larger J\f, the lower the relative dis- 
persion. In our case: 



<Yz) 
Yz 

CT(Miock) 



M, 



oc 



lock 



and N{Yz) cx 

and 7V(Afiock) oc Ms 



(15) 



and for the case of the true yield pz = Yz /Miock it can be 
demonstrated that (see Appendix in perviiio et al. 2002, 
for details): 



g(pz) 



and 7V(pz) oc Ms 



(16) 



Note that M has the advantage of scaling always with 
Ms whatever the considered quantities. 

With this formalism we have computed the evolution 
with time of the true yields and the corresponding M val- 
ues for Hydrogen, Helium, Carbon, Nitrogen and Oxygen 
at five metallicities (Z=0.020, 0.008, 0.004 and 0.0004). We 



have used the ejected masses presented in Portinari et al 
(1998) for massive stars (M> 8 M©). These yields include 
the contribution of stellar winds and those from the type II 
supernova explosions. For in termediate a nd low mass stars 
we have used the yields from Buell ( 1997 ), interpo lated for 
the above m etallicities, such as those shown in [Gavilan 



et al 



(|2002 | ) . We have considered as solar abundances 
(in mass) the values H=0.735, He=0.248, C= 2 .92 x lO"'^ 
N=8.56 X 10"'^ and 0=7.95 x 10"^ taken from ICrevesse & 



Sauvaj (1998). The abundances of Carbon, Nitrogen and 



Oxygen have been scaled linearly with the metallicity Z for 
the other metallicities. The Helium abundance has been 
considered as He=0.24 for metallicities different from so- 
lar. We had considered a Salpeter IMF slope (a — —2.35) 
in the mass range 120 to 0.08 M©. Metallicity dependent 
lifetimes have been taken from Portinari et al. ( |1998| ) for 
all the mass range. 

In Fig. [l] we show the evolution of the true yields 
of ph, PHc, PC: Po and and their relative dispersion 
as a function of time. The relative dispersion has been 
obtained including the covariance and the dispersion in 
A/lock- However, the values do not differ if it is consid- 
ered that Miock has no dispersion at all, except for ages 
lower than 10 Myr. The figure shows that N{po) reaches 
a constant value at ^ 50 Myr, when the last supernova 
explodes {mr{t = 50Myr) = 8 Mq). However, there is an 
evolution in the po value mainly due to the evolution of 
-A/lock, that decreases with time, and a small production 
of Oxygen by low mass stars (see below). The values of 
the true yields pc and pn and their associated J\f values 
increase with time. This reflects the contribution of low 
mass stars to the total yield, which result in an increment 
on the effective sources that produce these elements and, 
so, a lower relative dispersion (i.e. larger Af values). The 
effect is more significant in the case of Nitrogen, where 
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Z = 0.0004 
Z = 0.004 
Z = 0.008 
Z = 0.020 



Age [yr] Age [yr] 

Fig. 1. Evolution of the true yields pn, PHo, Pc, Po and pn and their relative dispersions as a function of time. 



the contribution of stars with a lifetime larger than 50 
Myr (initial masses lower than 8 Mq) produces a signifi- 
cant enhancement on the yield. However, the increment is 
similar for all the considered metallicities, contrary to the 
expectations for a secondary production of Nitrogen. On 
the other hand, the ejection rates used show an strong de- 
pendence of the Carbon yield with the initial metallicity, 
and hence showing a pseudo-secondary behaviour. 

We show the evolution of the correlation coefficient of 
p(yH,i%c), p(lH,ro), piYKe,Yo), piYo,Y^), p{Yo,Yc), 
and p{Yc,Y^) in Fig. ||. Note that the value of the covari- 
ance coefficient is an absolute one, independent of sam- 
pling effects. In this case it is a measure of the productions 
of the elements among different types of stars. 

The figure shows that 1h and Yq, and Yh and Yhc are 
almost anti-correlated for all the ages and metallicities. In 
a first order approximation, it means that the regression 
coefficient (i.e. the correlation coefficient) of a linear fit 
of the abundance of Helium or Oxygen vs. abundance of 
Hydrogen will be close to -1. In the case of p(1ho, ^o) the 
figure shows that a linear fit of the Helium abundance vs. 
12 + log (0/H), as the one used for the determination of 
the primordial Helium abundance, may have a value be- 
tween 0.9 and 0.7 (i.e. there will be an intrinsic dispersion 
in such type of representation of the data, as it is indeed 
observed). However, we want point out that we are work- 



ing with true yields, which are different from the observed 
ones. 

In the case of Yq and Fn the correlation is relatively 
strong (close to 1) in the first 50 Myr, when the produc- 
tion of Oxygen and Nitrogen is due to a similar population 
of stars. For ages larger than 50 Myr, the correlation co- 
efficient becomes dependent on the metallicity: the lower 
the metallicity, the closer to zero the correlation coefficient 
will be. A zero value in the correlation coefficient means 
that Oxygen and Nitrogen are produced by different types 
of stars. It also means, in a first order approximation (but 
see next section), that the N/0 ratio might have a larger 
dispersion in the low abundance regions, when the con- 
tamination of the ISM is due to stars with an initial low 
metallicity. 

The correlation between Carbon and Oxygen shows 
a different behavior. The correlation coefficient is also de- 
pendent on the initial metallicity, but in this case the lower 
value is 0.6, i.e. the production of Oxygen and Carbon is 
produced by similar (but not the same) stellar popula- 
tions. Also, the lower value of the correlation coefficient 
at solar metallicity implies that, in a first order approxi- 
mation, the dispersion in the C/0 ratio will increase with 
metallicity. 

Finally, we had also plotted the correlation between 
Carbon and Nitrogen, very similar to the one for Oxygen 
and Nitrogen. This is the expected behaviour if Carbon 
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Fig. 2. Evolution of the correlation coefficient p(yH, Ihc), p(1h, >b), piYsc, Yq), p(Yo, ^^n), p{Yo,Yc), and p{Yc,Yj^) 
as a function of time. Symbols like in Fig. |^. 



and Oxygen originate mostly from massive stars. Such a 
correlation coefficient was also shown in Cervifio et al. 
( |200l| ) for solar metallicity and ages between 1 and 10 
Myr both being quite similar in both cases despite the 
different treatment used. 

Once these true yields, their time evolution when a 
simple stellar population is created, and their dispersions 
have been computed, we can use them for calculating the 
abundances of these elements. As we have already ex- 
plained, we want to use simple models before treating 
this subject with numerical models, thus we will use these 
yields in the most simple GCE model: The Closed-Box 
Model. 



4. The Closed-Box Model 

The original formalism of the Closed-Box model was es- 
tablished by |Tinsley| ( |1980| ). pVIaeder) ( |1992| ) revised the 
formalism and we follow the prescriptions of this author. 
Let us assume that the total mass in the system, ill, is 
constant: M = s + g = c, where g is the amount of mass 
in the gas and s the amount of mass in stars. In this case 



Z{t) - Zo = pz In ■ 



-pz(t)lnM(0, (17) 



where fi is the gas fraction with an associated variance 
cr^ = crl/M'^, and = Z(i = 0) is the metalhcity at the 
onset of the burst of star formation. 

For simplicity, let us assume that there is only a disper- 
sion due to the discrete number of stars in Iz- In this case, 
the correlation coefficient p(lzi , ^Z2 ) equals to p(pzi , PZ2 ) • 

The variance in Z is then: 



(18) 



And the relative dispersion of the abundance ratio of 
two elements Zi respect to Z2 is: 



{Z1IZ2 



hi' /i 



zl 



'^p-z.^ 2cov(pzi,PZ2) 



Z\ Z2 



(19) 



Note that the dispersion scales with In /i. In a first 
order approximation it means that the dispersion on ratios 
of elements will be larger for low metallicity systems. 

The Closed-Box Model needs the use of an additional 
assumption about the true yield. In such model, the time 
evolution of the true yield is not considered, because, due 
to the hypotheses involved in its calculation, the true yield 
corresponds to the total yield obtained for a single stel- 
lar generation, that is, the maximum value obtained when 
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all stars die. This is the final value reached at the end of 
the time evolution in Fig. |l|. Thus, the abundance com- 
puted with the Closed-Box Model is an upper limit of this 
one which may be actually observed in a star forming re- 
gion. On other hand the Closed-Box model assumes that 
all stars are divided in two groups, a first one formed by 
the low mass stars living forever, and a second one with 
the more massive which will die immediately after being 
created. The mass limit between both groups depends on 
time. Therefore, to use adequately the Closed-Box model, 
it is necessary define a time scale (or a turn-off mass) r(m) 
where it is assumed that all the stars with lifetimes shorter 
than T(m) (or masses larger than m) produce all their 
yield instantaneously. Concerning the total mass trans- 
formed into stars, the Closed-Box model is not dependent 
on the assumed SFH, so, an amount of mass Mg trans- 
formed into stars or a constant star formation rate with a 
value Ms/T{m) produce the same result. Finally, it is also 
needed a value of ^ to describe the evolution of the abun- 
dance. In the following we will use a value of T{m) — 50 
Myr to obtain the dispersion in the age-0 /H relation and 
an value of T{m) ~ 1 Gyr to obtain the dispersion on the 
N/0 ratio. The metallicity evolution will be performed in 
terms of /i. Finally, we have assumed the same Zg values 
than the ones used for the computation of the true yields. 

4.1. Age-O/H relation 

Let us now consider the age-metallicity relation. Although 
this relation is usually given as a function of [Fe/H], we 



will use the abundance (0/H) as the metallicity indicator. 
To calculate the Iron yield it is necessary to include the 
ejection from type I SN explosions, which are not consid- 
ered in this work. 

The stellar population has produced most of the 
Oxygen at 50 Myr, so we have chosen this value of T(m) 
to obtain the dispersion in the age-0 /H relation. We have 
transformed the mass-abundance ratio in the element- 
abundance ratio by number obtaining this in the used way: 



(O/H) = 12- 



log — = 12 



log 



16 X Zh 



(20) 



The corresponding dispersion as a function of ^ is 
shown in the left panel Fig. |3[ The dispersion increases 
with lower values of fj, as is expected from Eq. (|l^) . 

As an example, let us assume the yield at Z=0.0004. 
It reaches a value of l2 + log{0/H) = 8.93 at ^ = 0.4. Let 
us also assume a constant star formation rate of 'i'(t) =4 
Mq pc~^ Gyr~^ for the solar cylinder (c.f. Pagej 1977, 
Table 7.7). In this situation, a dispersion oia^o/H) =28.5 
pc~^ dex is expected, i.e. ~ 0.3 dex if the star formation 
takes place in regions of 100 x 100 pc^. 

This value is similar to the observational errors. In gen- 
eral, the dispersion will scale with the amount of gas trans- 
formed into stars or with the star formation rate (once 
defined a time scale). So, in this simple framework, the 
scatter in the age-metallicity relation is related with the 
Star Formation History itself. Unfortunately, this state- 
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ment can be only verified with more sophisticated GCE 
but that is not the goal of this paper. 

It is also interesting to note that the recompila- 
tion of correlation coefficients of the fit of O/H versus 
Galactocentric distance for the MWG presented in Henry 
& Wo: they (199£) has a composite value of -0.63. In 
the case of the true yields, the correlation coefficient be- 



Finally, we want to note that in the case of GCE mod- 
els attempting to trace the mctallicity distribution of our 
Galaxy, the dispersion will not scale with the total mass 
of the Galaxy. Instead, it must be computed with the 
masses of the subsystems that are sampled: the disper- 
sion of abundances for individual Hii regions from the 



tween and Yq varies between 0.75 (Z^Z©) and 0.95 
(Z=0.004) with a strong dependence on metallicity. In our 
framework is not possible to compute the correlation coef- 
ficient of the relation between O/H and the Galactocentric 
distance, but a more elaborated model must be able to 
obtain a theoretical value of this correlation coefficient 
comparable with the observed one. 



4.2. N/0 ratio 

For the N/0 ratio, we had chosen a value of T(m) = 1 
Gyr that implies that all the stars with mass larger than 
~ 2 M© have enriched the ISM with their Nitrogen. The 
N/0 ratio has been computed with the relation: 



global metallicity gradient observed in our Galaxy ( Henry 
Worthey 1999| ) could be understood in this framework 



{N/0) 



1 ^ 



log 



14 X Z, 



o 



16 X Z 



N 



(21) 



The corresponding dispersion as a function of /i is 
shown in the right panel of Fig. |^. Again, a direct com- 
parison with observed data can not be performed. First 
of all, the {N/0) value is fixed by the yields used (and 
the T(m) value). However, it can be shown that a <J(n/o) 
value larger than 0.5 dex will be a natural effect if the 
ISM has been contaminated by a burst of star formation 
with an amount of gas transformed into stars lower than 
8 X 10^ M© or, equivalently, star formation rates lower 
than -^{t) = 80 M© Myr^^ (using the yield for Z=0.0004). 
Again, a more sophisticated model is needed to establish 
any conclusion. 

5. Conclusions 

In this work we have shown a theoretical formalism for 
the evaluation of the dispersion on GCE models. Using 
a simple Closed-Box model, we have obtained a first or- 
der approximation to the relevance of the discreteness of 
the stellar population and the sampling effects in GCE 
models. Despite the approximations, the effect of the dis- 
creteness of the stellar populations may be relevant for 
explaining the dispersion in the age-metallicity relation. 
Such observed dispersion may be an indicator of the SFH 
of our Galaxy. It may be also relevant in the study of the 
observed dispersion in the N/0 ratio in Hii regions and 
dwarf galaxies. 

We have found that the dispersion in the N/0 ratio is 
larger for lower metallicities. This effect is indeed observed 
in our Galaxy. We have obtained a first order theoretical 
estimation for the goodness of a linear fit of the Helium 
abundance vs. 12 -I- log (O/H) with values of the regression 
coefficient between 0.9 and 0.7. 



since these regions are less massive systems. 

However, no conclusion can be obtained unless this ef- 
fect is included in a more sophisticated GCE. It is there- 
fore necessary to make a more exhaustive study on this 
subject. It will be performed in forthcoming papers. 
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